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We employ non— perturbative flow equations to compute the equation of state for 
two flavor QCD within an effective quark meson model. Our treatment covers 
both the chiral perturbation theory domain of validity and the domain of validity 
of universality associated with critical phenomena. In the vicinity of T c and zero 
quark mass we obtain a precision estimate of the universal critical equation of 
state of the three dimensional 0(4) symmetric Heisenberg model. For realistic 
quark masses the pion correlation length near T c turns out to be smaller than its 
zero temperature value. 



QCD in a thermal equilibrium situation at sufficiently high temperature 
differs in important aspects from the corresponding zero temperature or vac- 
uum properties u. A phase transition at some critical temperature T c or a 
relatively sharp crossover may separate the high and low temperature physics. 
Concentrating on the chiral aspects of QCD the transition is .u&Lated to a qual- 
itative change in the chiral condensate. It was pointed out cm that for suffi- 
ciently small up and down quark masses, m u and mj, and for a sufficiently 
large mass of the strange quark, m s , the chiral transition is expected to belong 
to the universality class of the three dimensional 0(4) Heisenberg model. It 
was suggestedu that a large correlation length max lead to a disoriented chiral 
condensate u with possible distinctive signatures 13 in a relativistic heavy-ion 
collision. The question how small m u and rrid would have to be in order to see 
a large correlation length near T c and if this scenario could be realized for real- 
istic values of the current quark masses remained, however, unanswered. The 
reason was the missing link between the universal behavior near T c and zero 
current quark mass on one hand and the known physical properties at T = 
for realistic quark masses on the other hand. Lattice QCD seems particularly 
suitable for such a study, however, exploring the universal region is limited by 
present computer resources El. 

It is the purpose of this talk to provides the "missing link" . Our approach 
is based on the use of a non-perturbative flow equation for a scale dependent 
effective action which is the generating functional of the 1PI Green 

"Talk given at Eotvos Conference in Science: Strong and Elcctroweak Matter, Eger, 
Hungary, 21-25 May 1997. 
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functions in the presence of an infrared cutoff ~ k. Varying the infrared cutoff 
k allows us to consider the relevant physics in dependence on some momentum - 
like scale. The standard effective action is obtained by removing the infrared 
cutoff (k — > 0) in the end. The fcrdependence of the effective average action is 
given by an exact flow equation 0, which for scalar fields <&i reads 

(2) 

Here 1^ denotes the matrix of second functional derivatives of Tj, with re- 
spect to the field components. We use a momentum dependent infrared cutoff 
Rk(q) = Z$>^q 2 e~ q ~ l k /(l — e~ q l k ) with an appropriate wave function 

renormalization constant. In momentum space the trace reads Tr = J Y]j- 
We employ for scales below a "compositeness scale" of fc$ ~ 600 MeV a 
description in terms of quark and scalar mesonic degrees of freedom. This effec- 
tive quark meson model can be obtained from QCD in principle by "integrating 
out" the gluon degrees of freedom and by introducing fields for composite op- 
erators 0. In this picture the scale fc$ is associated to the scale at which the 
formation of mesonic bound states can be observed B in the flow of the mo- 
mentum dependent four-quark interaction. We imagine that all other degrees 
of freedom besides the quarks ip and the scalar and pseudoscalar mesons con- 
tained in the complex field $ are integrated flut. Our truncation corresponds 
to the ansatz for the effective average actionQ 

T k = J d 4 xl Z^ a i^ a + Z $ , fc tr [d^d"*] + lT fc (4,*t) 

(2) 

+W (^*-«- " ^(* f u) ip" - ^ tr ($tj + /$) 

Here is invariant under the chiral flavor symmetry SUl(2) x SUr(2) in 
absence of the explicit symmetry breaking through the source term j ~ m = 
(rn u + rrid)/2. We will consider the flow of the most general form of the po- 
tential term Uk consistent with the symmetries. For k — > the potential Uk 
encodes the equation of state. At non-zero temperature the ansatzB therefore 
allows to study the complete non-analytic behavior of the effective potential, 
or equivalently the free energy, near the critical temperature of the second or- 
der phase transition. On the other hand, our approximations for the kinetic 
terms are rather crude and parameterized by only two running wave function 
renormalization constants, Z$k and Z^±. The same holds for the effective 
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Yukawa coupling h^. We further neglect the scalar triplet clq and the pseu- 
doscalar singlet (associated with the 77') for k < fc$ ||. This can be achieved in 
a chirally invariant way and leads to the 0(4) symmetric linear sigma model 
for the pions and the sigma resonance, however, coupled to quarks now. The 
non-perturbative flow equations for the quark-meson model are obtained from 
eq. ([!]) generalized □ to include fermions and using the ansatz (||) for Tp. . We 
solve them numerically. 

A reliable quantitative derivation of the effective quark meson model from 
QCD is still missing. We emphasize, however, that the quantitative aspects 
of this derivation will be of minor relevance for our practical calculations in 
the mesonic sector: If the effective Yukawa coupling between the quarks and 
the mesons turns out to be strong at the compositeness scale we observe a 
fast approach of the scale dependent effective couplings to approximate partial 
infrared fixed points As a consequence, the detailed form of the meson 
potential at becomes unimportant, except for the value of one relevant scalar 
mass parameter . Here we fix m*^ from phenomenological input such that 
f„ = 92.4 MeV (for m n — 135 MeV) which sets our unit of mass for two flavor 
QCD. The only other input parameter we use is the constituent quark mass M q 
to determine the scale k$. We consider a range 300 MeV 55 M q <, 350 MeV and 
find a rather weak dependence of our results on the precise value of M q . We 
point out that though a strong Yukawa coupling at fc$ is phenomenologically 
suggested by the comparably large value of the constituent quark mass M q it 
enters our description as a (consistent) assumption. 

The equation of state expresses (ipip) as a function of TJ-ajid rh where the 
chiral condensate is related to the expectation value ($) byBfl 

(W) = -2m 2 k J($)-m], (3) 

Curve (a) of figure |l| gives the temperature dependence of (f/'V') m the chiral 
limit rh = 0. Jfcijc the lower curve is the full result for arbitrary T whereas 
the upper curve&til corresponds to the universal scaling form of the equation of 
state for the 0(4) Heisenberg model. We see perfect agreement of both curves 
for T sufficiently close to T c = 100.7 MeV. This demonstrates the capability 
of our method to cover the critical behavior and, in particular, to reproduced 
the critical exponents of the 0(4)-model. The curves (6), (c) and (d) are 
for non-vanishing values of the average current quark mass rh. Curve (c) 
corresponds to m p h ys or, equivalently, m^T = 0) = 135 MeV. One observes 
a crossover in the range T — (1.2 — 1.5)T"c- In order to facilitate comparison 

6 The present investigation for the two flavor case does not take into acepunt a speculative 
"effective restoration" of the axial !7a(1) symmetry at high temperature □. 
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Figure 1: The plot shows the chiral condensate ^V 1 /") as a function of temperature T. Lines 
(a), (6), (c), (d) correspond at zero temperature to = 0, 45 MeV, 135 MeV, 230 MeV, 
respectively. For each pair of curves the lower one represents the full T— dependence of (ipip) 
whereas the upper one shows for comparison the universal scaling form of the equation 
of state for the 0(4) Heisenberg model. The critical temperature for zero quark mass is 
T c = 100.7 MeV. The chiral condensate is normalized at a scale ~ 620 MeV. 



with lattice simulations which are typically performed for larger values of 
we also present results for (T = 0) = 230 MeV in curve (d) . One may define 
a "pseudocritical temperature" T pc associated to the smooth crossover as the 
inflection point of (tpip) (T). The value for the pseudocritical temperature for 
m„ = 230 MeV is T pc ~ 150 MeV. For realistic quark mass, or m v — 135 MeV, 



we obtain T, 



pc 



L pc 

130 MeV. 



A second important result of our investigations is the temperature depen- 
dence of the space-like pion correlation length to^T^T). The plot for m w (T) 
in figure ^| again shows the second order phase transition in the chiral limit 
rh = 0. In this limit the behavior for small positive T — T c is characterized by 
the critical exponent v, i.e. m w (T) = (£ + ) T c ((T — T c )/T c ) and we obtain 
v = 0.787, £ + = 0.270. For rh > we find that m,(T) remains almost constant 
for T Si T c . For T > T c the correlation length decreases rapidly and for T>T C 
the precise value of rh becomes irrelevant. The overall size of the pion correla- 
tion length near the critical temperature is given by m^(T pc ) ~ 1.7m ff (0) for 
the realistic value m p h ys - 

We point out two important answers one obtains from this study: First of 
all, for a thermal equilibrium situation the chiral transition gives no indication 
for strong fluctuations of pions with long wavelength fl. The longest correla- 
tion lengths near the crossover temperature T pc even becomes smaller than at 
T = On The second answer concerns the applicability of universal, i.e. almost 

c See ref. i for a discussion in the context of a chiral transition far from equilibrium. 
d li should be emphasized, however, that a tricritical behavior with a massless excitation 
remains possible for three flavorslij. 
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Figure 2: The plot shows m n as a function of temperature T for three different values of 
the average light current quark mass rh. The solid line corresponds to the realistic value 
rh = rhphys whereas the dotted line represents the situation without explicit chiral symmetry 
breaking, i.e., m = 0. The intermediate, dashed line assumes rh = m p i lys /10. 



model independent, arguments for a description of the chiral phase transition. 
Despite the observed comparably short correlation lengths at non-zero tem- 
perature, we find the approximate validity of the 0(4) scaling behavior over a 
large temperature interval near and above T c even for quark masses somewhat 
larger than the realistic ones. The_present approach can be extended to the 
case with three light quark flavorsll^. 
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